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Abstract 

We examine elastic Compton scattering from the deuteron for photon energies 
ranging from zero to 100 MeV, using state-of-the-art deuteron wave functions and 
AA-potentials. Nucleon-nucleon rescattering between emission and absorption of the 
two photons is treated by Green’s functions in order to ensure gauge invariance and the 
correct Thomson limit. With this Green’s-function hybrid approach, we fulfill the low- 
energy theorem of deuteron Compton scattering and there is no significant dependence 
on the deuteron wave function used. Concerning the nucleon structure, we use Chiral 
Effective Field Theory with explicit A(1232) degrees of freedom within the Small 
Scale Expansion up to leading-one-loop order. Agreement with available data is good 
at all energies. Our 2-parameter fit to all elastic yd data leads to values for the static 
isoscalar dipole polarizabilities which are in excellent agreement with the isoscalar 
Baldin sum rule. Taking this value as additional input, we find = (11.3T0.7 (stat)± 
0.6 (Baldin))-10“^ fm^ and Sm — (3.2=f 0.7 (stat)±0.6 (Baldin))-10“^ fm^ and conclude 
by comparison to the proton numbers that neutron and proton polarizabilities are 
essentially the same. 
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1 Introduction 


Analyzing protons and neutrons with electromagnetic probes has a long history in the held of 
nucleon-structure studies. In Compton scattering off a single nucleon, the electromagnetic 
held of the scattered photon attempts to deform the nucleon. The global resistance of 
the nucleon’s internal degrees of freedom against this deformation is measured in terms 
of the electromagnetic polarizabilities, which makes them an excellent tool to study the 
structure of the nucleon. In principle, polarizability is a function of the frequency of the 
electromagnetic wave. Therefore, energy-dependent polarizabilities have been introduced in 
Refs. [21 El- In this work, we determine the static values of the nucleon polarizabilities from 
experiment, i.e. the values in the limit of vanishing photon energy, which we therefore denote 
as the polarizabilities for simplicity. The two most prominent nucleon polarizabilities are the 
static electric and magnetic dipole polarizabilities and Pm- For the proton, rather precise 
experimental values for these two quantities exist, e.g. those given in a recent review | 1 ], 
which were obtained as the weighted average over several experiments, 

= (12.0 ± 0.6) ■ 10"^ fm^ = (1.9 ^ 0.6) ■ 10“^ fml (1.1) 

Our recent, Baldin-constrained £t to the proton Compton data yielded [ 3 ] 

= (11.04 ± 1.36(stat) ± 0.4(Baldin)) ■ lO”'^ fm^, 

= (2.76 =F 1.36(stat) ± 0.4(Baldin)) ■ lO”'^ fm^, (1.2) 

in agreement with Eq. (dD The framework that has been chosen for this extraction is 
third-order Small Scale Expansion (SSE), which is one possible way to include the 

A(1232) resonance explicitly in Chiral Effective Field Theory, cf. Ref. [B]. It is this very 
framework on which the one-nucleon sector of our present work is built. 

On the other hand, neutron polarizabilities are much harder to access experimentally, 
as there is no stable single-neutron target for Compton scattering. Therefore, one has to 
rely on other methods in order to investigate these quantities. One approach is quasi-free 
Compton scattering off the neutron bound in the deuteron, or the scattering of neutrons 
from a lead target. The weighted average over several experiments investigating these two 
processes gives the result ji] 

= (12.5 ± 1.7) ■ 10"^ fm^ = (2.7 ^ 1.8) ■ 10“^ fml (1.3) 

The numbers given in Eq. (Q do not include information from elastic low-energy Comp¬ 
ton scattering from light nuclei, as these processes so far led to values inconsistent with 
Eq. fll.3j) . However, several such experiments have already been performed jZl El El and 
further proposals exist, e.g. from the deuteron or ^e at TUNL/HI 7 S ^U] and on deuteron 
targets at MAXlab CH as well as at the S-DALINAC ^21- They promise an extensive study 
of elastic deuteron Compton scattering below the pion-production threshold. From a theo¬ 
rist’s point of view, extracting the neutron polarizabilities requires an accurate description 

^The central value of the Baldin sum rule + ~ (13.8±0.4) • 10“^ fm^ has been used in Eq. IE21. 


2 



both of the structure of the nucleon and of the dynamics of the low-energy degrees of free¬ 
dom within the deuteron, as one has to account for the nucleon polarizabilities as well as for 
meson-exchange currents. We remind the reader that the deuteron as the proton-neutron 
bound state is an isoscalar object, so that only a% = ^ {a^ + a'%), | + Plea¬ 
sured. The isovector polarizabilities are dehned as = (ozg —a^), | . A 

hrst attempt to £t the isoscalar polarizabilities to the elastic deuteron Compton-scattering 
data from Illinois [7] and SAL jH] was made in The extracted neutron polarizabilities 

= (9.0 ± 3.0) ■ 10“^ fm^, [3^ = (11.0 ± 3.0) ■ lO”'^ fm^ indicated the possibility of a 
rather large isovector part, especially in the magnetic dipole polarizability, in contrast to 
the quasi-elastic result from m contained in Eq. (jl.3|) . The enhancement of j3^ in jT^j is in 
our opinion due to an insufficient description of the SAL-data around 94 MeV, in particular 
in the backward direction, which have been a puzzle for several years. It was hnally resolved 
by the authors of P] in an Effective-Field-Theory calculation by observing that dynamical 
effects from explicit A(1232) resonance degrees of freedom are large at these energies. The 
calculation presented in P agrees reasonably with data but only works above some lower 
energy limit of the order 50-60 MeV. The “best” (Baldin-constrained) fit results for the 
isoscalar polarizabilities given in pp are in excellent agreement with Eqs. dHD and (d; 

a% = (12.6 ± 0.8(stat) ± 0.7(wf) ± 0.6(Baldin)) ■ 10“^ fm^, 

PIj = (1.9 =F 0.8(stat) =F 0.7(wf) ± 0.6(Baldin)) ■ 10“^ fm^, (1.4) 

with ’wf’ denoting the uncertainty arising from the wave-function dependence observed 

in p. 

On the theory side. Heavy Baryon Chiral Perturbation Theory (HByPT) predicts that 
the proton and neutron polarizabilities are equal at leading-one-loop order ra- Comparing 
Eq. (jl.lj) to ()1.3j) or ()1.4j) . respectively, suggests that as and {3m have indeed only small 
isovectorial components. 

In the present work, we investigate elastic deuteron Compton scattering including photon 
couplings to the one-pion exchange. Diagrams which are characterized by the propagation 
of the nucleons between the two photon interactions are calculated using Green’s-function 
methods in analogy to Ref. unj. The single-nucleon structure is included within Chiral 
Effective Field Theory (y;EFT) with explicit A degrees of freedom up to leading-one-loop 
order like in Ref. P, where isovectorial components are zero. The A(1232) is treated 
according to the power-counting rules of the Small Scale Expansion p. We aim for a 
consistent calculation at photon energies u below 100 MeV, i.e. we 

1 ) give an improved description of the high-energy (m 95 MeV) data from P with 
respect to the calculations presented in naiini. 

2) present an alternative approach that does not suffer from the shortcomings of the EFT 
calculations of Refs. pimiiHi, which are not applicable below some lower energy limit. 
Stated differently, we demand that our calculation reaches the correct static limit, i.e. 
the Thomson limit for deuteron Compton scattering. 
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3) contribute to the ongoing discussion of the neutron polarizabilities via hts of the 
isoscalar polarizabilities a% and to all existing elastic deuteron Compton data. 

Effective Field Theory has already been used for the latter purpose in pp and ^H]- The 
authors of iniiHi, as well as those of Ref. na. followed Weinberg’s proposal na to cal¬ 
culate the irreducible kernel for the 'yNN —>■ yiViV process in Effective Field Theory and 
then folded this with external deuteron wave functions, derived from high-precision NN- 
potentials such as Nijm93 [201, CD-Bonn |2I] or AV18 j23- This “hybrid” approach has 
proven quite successful in describing scattering reactions like nd G~d |21] and other 
processes, see e.g. also |23]. However, the C>(p^)-HBxPT calculation of Ref. [TH] in the 
Weinberg hybrid approach gives an insufficient description of the SAL-data [2] measured 
in the backward direction. Therefore, the authors of excluded these (two) data points 
in some of their fits. Those of Ref. PP had to restrict themselves to the published data 
above 60 MeV. It is one of the central results of the approach presented in this work that 
we do not need any such constraints, in particular in view of forthcoming data even beyond 
w = 100 MeV dH. 

It is worth noting that a fully self-consistent xEFT of the two-nucleon system is not yet 
available. To take the counting rules of the one-nucleon system simply over to the potential 
of the two-nucleon system was Weinberg’s original proposal m- However, Nogga et ah and 
Birse dOEII showed recently that this approach is not self-consistent in W-scattering, 
but that additional short-distance interactions are needed in higher partial waves to recover 
phase shifts which are insensitive to details of short-distance physics. As this problem does 
not affect the ^Si and ^Di channels, Weinberg’s suggestion [12] of a hybrid approach seems 
still to be applicable in deuteron Compton scattering. 

In this work, we present a new hybrid approach, including the full two-nucleon Green’s 
function in all diagrams with an W-intermediate state, which was only treated perturba- 
tively in Ref. PP, according to the power-counting rules of third-order SSE for high-energy 
external probes. Besides the single-nucleon current, we couple the photon held also to the 
meson-exchange currents of the two-nucleon system. The strictly perturbative calculations 
of [TJ [T7| dHI have the disadvantage that they become inapplicable below oj ~ 50-60 MeV. 
They even over-predict the deuteron Thomson limit by more than a factor 2 |2H|- This is 
the more damaging as the Thomson limit is a simple consequence of gauge invariance, so 
that the calculation in such a simplistic application of Weinberg’s counting rules obviously 
violates gauge invariance. The calculation presented in this publication does not suffer from 
that limitation. In fact, it reaches the exact Thomson limit. Consequently we are able to 
describe the low-energy (a; < 60 MeV) data well. The various cross-checks we perform in 
this presentation make us conhdent that the cross sections and polarizabilities we obtain in 
this paper will be essentially unchanged once more systematic theories of nuclear interac¬ 
tions are found. Its relation to a more rigorous EFT approach will be elaborated on in a 
future publication on power-counting in two-nucleon xEFT Q. 

Comparing to Refs. usiiinj, we see the main advantage of our calculation in an improved 
treatment of the single-nucleon Compton multipoles, which in those works are only included 
via the leading uni and subleading ra terms of a Taylor expansion in oj, whereas we keep the 
full energy dependence of the Compton multipoles, including the explicit A(1232) following 
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third-order SSE as worked out in The huge influence of the A(1232) in single-nucleon 
Compton backscattering is well-known. It is due to the strong paramagnetic Ml coupling 
of the photon to the N —>■ A transition. A similarly strong influence of the A resonance has 
been found in deuteron Compton scattering in the backward direction in pQ. Therefore, we 
advocate to retain this degree of freedom explicitly in any Compton calculation. 

This paper is organized as follows: In Section we give a brief survey of the theoretical 
formalism applied, focussing on those diagrams which have a two-nucleon intermediate state. 
We close the section by demonstrating that our approach is gauge invariant and therefore 
fulfills the low-energy theorem EDI. i.e. it generates the correct static limit. In Section |2l we 
discuss our predictions for four different photon energies, ranging from 50 MeV to 100 MeV, 
and we compare those to data and to the 0{e^) SSE calculation of Ref. [Tj. The subsequent 
Section 0] contains our fits of the isoscalar polarizabilities to all existing elastic deuteron 
Compton data. We show that the data are self-consistent and in good agreement with the 
theoretical expectation that isovector components are small. We conclude in Sectional hav¬ 
ing shifted the most technical parts to Appendices El and ^ and to Ref. ESI. In Appendix O 
we derive the total deuteron-photodisintegration cross section from our Compton amplitude 
in order to cross-check our calculation. We And perfect agreement with data and predictions 
from Effective Range Theory. 

2 Theoretical Framework 

In this work, elastic deuteron Compton scattering is examined, including explicit A reso¬ 
nance degrees of freedom. In general, the T-matrix for Compton scattering off the deuteron 
is derived as the matrix element of the interaction kernel, evaluated between initial- and 
final-state deuteron wave functions, as explained in great detail in in|: 

T = (I'/l/C, + A', a j) (2.1) 

The first piece in Eq. m is called the “two-nucleon irreducible” part and the second is 
the “two-nucleon reducible” part. Two-nucleon reducible diagrams are those which contain 
an intermediate state with only the two nucleons as particle content. G is the two-particle 
Green’s function, constructed from the two-nucleon irreducible interaction V and the free 
two-nucleon Green’s function. denotes the coupling of one photon to the two-nucleon 
system, is the two-nucleon irreducible kernel for the coupling of the two photons. 

The main difference between this work and Ref. pQ is the treatment of those diagrams 
which are characterized by the propagation of the two nucleons in the intermediate state 
between the two photon interactions. In nnmiiHi, such diagrams have been calculated 
following the power-counting rules of Effective Field Theory for “high-energy” external 
probes, i.e. u ~ m^. For large photon energies, a perturbative treatment is possible as is 
easily understood heuristically: The absorption of a high-energy photon immediately tears 
the two nucleons apart, so the deuteron would be destroyed if the second photon was not 
emitted near-instantaneously. It turns out that up to leading-one-loop order in HByPT, as 
well as in SSE, the only diagrams with a two-nucleon intermediate state are the nucleon- 
pole diagrams, sketched in Fig. E Note that the s-channel diagram, i.e. the left diagram 
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in Fig. Q is part of K-yG K^, as it has an intermediate state with only the two nucleons 
propagating. While the nucleons can in general interact with each other between the two 
photon vertices, the authors of Refs. mini calculated only up to leading-one-loop order, 
where no such interactions are possible. Therefore, we sketch the propagation of two free 
nucleons in Fig. ^ 



Figure 1: Nucleon-pole terms without rescattering of the nucleons between the two photon 
interactions. The photon-nucleon vertex is given by minimal substitution or magnetic- 
moment interactions. 


As a consequence of the power counting applied, the calculations of Refs. BimiiHi 
break down in the low-energy regime, manifested in an incorrect Thomson limit. Their lower 
energy limit was found to be of the order of 50-60 MeV. The reason for the mismatch at low 
energies is that one has to treat the np-interaction non-perturbatively. Therefore, we now 
include the full np-S'-matrix in the intermediate state. The possible rescattering between 
the two nucleons is denoted by a square in Fig. |21 where we sketch the differences between 
Ref. PP and this work. Foremost, we do not only include the free propagation of the two 
nucleons in the nucleon-pole diagrams. Fig. Q but we construct the full np-Green’s function 
whenever a two-nucleon intermediate state is involved. Furthermore, we allow for more 
ways of coupling the photon held to the two nucleons with respect to Refs. mini: Besides 
coupling to the single-nucleon current, like in Fig.^ we also include pion-exchange currents 
as shown in the lower line of Fig.|21 making use of Siegert’s theorem |2Ij, cf. Section lT^ Such 
diagrams appear only at 0{p'^) in HByPT and are therefore not included in pp. In Ref. |18j . 
the np-rescattering in the intermediate state is not fully included. It is well-known that 
only full inclusion of np-rescattering, together with the appropriate explicit pion-exchange 
diagrams, generates the Thomson limit of deuteron Compton scattering [13113 Ea ESI as 
direct consequence of only demanding gauge invariance |d()j . 

In this publication, we follow closely the work of Karakowski and Miller to construct 
the rescattering and the photon coupling to the meson-exchange currents via Siegert’s the¬ 
orem, cf. Fig. la However, before we turn to the calculation of the diagrams with an 
intermediate np-state, we recall all other contributions, which therefore are part of 
see Eq. (jnj. As those terms have already been discussed in Ref. PP and partly in earlier 
references therein, we shall be brief in the upcoming section. 
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Figure 2: Upper line: Sketch of the different treatment of diagrams with two-nucleon in¬ 
termediate state in the present work with respect to Refs. dim The square symbolises 
the NN-S-matiix. The blobs denote photon coupling to the one-body current and possible 
one-pion exchange, as indicated in the lower line. 


2.1 Diagrams without Intermediate Scattering 

Except for the diagrams with two-nucleon intermediate states, see Fig. El we apply the 
power-counting rules of the Small Scale Expansion, where the expansion parameter is e, 
denoting either a small momentum, the pion mass or the mass difference Aq between the 
real part of the A mass and the nucleon mass. We refer the interested reader to 133 
for the Ntt Lagrangean and to |SS! and jHj for the relevant pieces of the Avr Lagrangean. 
The power-counting scheme that we use for the nucleon-structure part of our calculation is 
motivated by Weinberg’s idea to count powers only in the interaction kernel ng. This hybrid 
approach is a well-established tool by now. While the kernel is power counted according to 
the rules of the Effective Field Theory, a deuteron wave function from state-of-the-art NN- 
potentials is used. Unfortunately, a fully self-consistent yEFT of the two-nucleon system 
is not yet available |2nil2Zl- We apply therefore here another hybrid approach which uses 
basic facts of nuclear phenomenology in the two-nucleon sector such as np-rescattering and 
meson-exchange currents. Its relation to a more rigorous EFT approach will be given in a 
forthcoming work on power counting in two-nucleon x^FT [21]. Still, the deuteron wave 
functions that we use are derived from modern AA-potentials: the AV18-potential ra and 
the NNLO chiral potential (SHI; where this last potential also follows Weinberg’s suggestion 
and is derived by applying HByPT power counting to the iViV-potential V. 

Considering only contributions without intermediate two-nucleon states, there is another 
possibility to classify diagrams. It is the separation into one-body and two-body pieces, 
i.e. into diagrams where only one of the two nucleons or both of them are involved in 
the Compton-scattering process. Obviously, such a rigorous separation is no longer possible 
when we calculate the diagrams including iViV-rescattering between the photon interactions. 
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cf. Fig. 121 As in Ref. P], we calculate up to 0{e^), i.e. there is no difference in the 
diagrams without an intermediate np-state between this work and pp. Therefore, we refer 
to Ref. PP for further details on these contributions and only list again the various diagrams 
for completeness: 

• One-body contributions without explicit A(1232) degrees of freedom. These are the 
single-nucleon seagull, Fig. Ha), and the contributions from the leading chiral dynam¬ 
ics of the pion cloud around the nucleon (Figs. Hc)-(f))- The pion pole (Fig. Hb)), 
i.e. the Tr'^-exchange in the t-channel, does not contribute to deuteron Compton scat¬ 
tering at this order, as it is isovectorial. Note that the nucleon-pole terms. Fig. ^ also 
contribute at 0{e^) PP. In this work, however, these diagrams are not only included 
perturbatively, like in Ref. pP , but we calculate them using the full np-Green’s function 
G, cf. Eq. fl2.1|l . Their evaluation is postponed to Section YI7]\ 

d 

(a) (b) (0) (d) (e) (f) 

Figure 3: One-body interactions without A(1232) propagator contributing to deuteron 
Compton scattering up to 0{e^) in SSE. Permutations and crossed graphs are not shown. 



• One-body diagrams with explicit A degrees of freedom, as shown in Fig.^J The A-pole 

diagrams (Fig. H^)) ^^e contributions from the pion cloud around the A(1232) 

(Figs.|li;b)-(e)). 

• Two isoscalar short-distance one-body operators (Fig. HKf)), whose coupling constants 
are denoted as pii 7 , png in Table [TJ These operators, which we determine via hts to ei¬ 
ther proton or deuteron Compton cross-section data, contribute energy-independently 
to the dipole polarizabilities a% and They are formally of 0{e^) but turn out 
to give an anomalously large contribution to the single-nucleon Compton amplitude. 
Therefore, we promote them to next-to-leading order as discussed in detail in [3]. 

• Two-body contributions with one pion exchanged between the two nucleons (Fig. E)). 
As discussed in Ref. P], the meson-exchange diagrams are identical in third-order 
HB^PT and SSE. These diagrams, together with those given in Figs. |2l and cire 
responsible for complying with the low-energy theorem, as discussed in Section 12.2.41 


All these diagrams (Figs. OHSI) are included in our interaction kernel. The SSE single¬ 
nucleon amplitudes can be found in [3], while the two-body contributions are given explicitly 
in HU- The only diagrams that remain to be calculated are the diagrams with an interme¬ 
diate np-state, cf. Fig. El 
















(b) 


(c) 


(d) 


(e) 


(f) 






Figure 4: Additional one-body interactions which contribute to deuteron Compton scatter¬ 
ing at 0{e^) in SSE compared to third-order HByPT. Permutations and crossed graphs are 
not shown. 



(a) (b) (c) (d) (e) 


Figure 5: Two-body interactions contributing to the kernel for deuteron Compton scattering 
at 0{e^) in SSE. Diagrams which differ only by nucleon interchange are not shown. 


2.2 Diagrams with Intermediate A/^A/^-Scattering 

2.2.1 A Short Note on the Thomson Limit 


Before we describe how we calculate the AA-rescattering diagrams sketched in Fig. |2l we 
consider the well-known Thomson limit, i.e. the deuteron Compton amplitude in the limit 
of vanishing photon energy: 


A 


Thomson 

d 






( 2 . 2 ) 


Friar showed that Eq. (O is a consequence of gauge invariance EDI. i.e. reaching the correct 
static limit is equivalent to gauge invariance as long as there are no photons in internal 
loops. Therefore, the EFT calculations miiTiiiHi, which are by construction inapplicable 
in that limit, would there also violate gauge invariance, albeit their interaction kernel is 
gauge invariant by construction. The violation appears when evaluating the kernel between 
the deuteron wave functions, without allowing the two nucleons in the intermediate state 
to interact with each other. The reason is that the deuteron wave function implies this 
interaction, which can be interpreted as the exchange of mesons, e.g. of pions, between the 
two nucleons. In order to achieve gauge invariance, it is therefore mandatory to include 
rescattering of the two nucleons on one hand and to couple the photons to these meson- 
exchange currents, cf. Fig. Inland Refs. [TIH ITT)] . It is one of the main advantages of this 
work with respect to [niniiiHi that our calculation does fulhll Eq. (El. 
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Reaching the Thomson limit is a non-trivial check because the deuteron mass is involved, 
whereas the Thomson seagull for Compton scattering from the proton, Fig. Ha), yields the 
amplitude 


a: 


Thomson 




rur. 


(2.3) 


with rUp the proton mass. The single-neutron amplitude vanishes in the static limit. There¬ 
fore, all other contributions to deuteron Compton scattering in the limit a; —*■ 0 have to 
cancel half of the proton amplitude (jZSD. In Section 12.2.41 we discuss how this cancellation 
comes about, after sketching in Sections 12 .2.21 and l2.2.dl onr way to calculate diagrams with 
intermediate np-scattering. 


2.2.2 Dominant Terms 


In the following, we briefly explain the “Green’s-function method” to include the NN- 
rescattering in the diagrams given in Fig. |2l using second-order time-ordered perturbation 
theory in the two-photon interaction. For further details see Appendix |B] and Refs. [THll28j . 
In general, the scattering amplitude for these processes can be written as 


c 


f ( df^ 7/ 

^int 

C){C 

Hint 


1 


-Ec 



( df, 7/ I I C, 7/, 7i) (C', 7/, % I I dj, % ) 

-^ + £;-b-Ec-^ 


(2.4) 


where we sum over all possible intermediate eigenstates C of the A^A^-interaction Hamilto¬ 
nian. The two terms in Eq. (El correspond to the two diagrams shown on the right hand 
side of Fig. 1^ We sketch these diagrams once again in Fig. IHlin order to explain the various 
terms in the denominators of Eq. (El, constituting the energy of the intermediate nucleons. 
In the yd cm frame, the incoming and outgoing photons have the same energy oj. B denotes 
the deuteron binding energy, is the kinetic energy of the incoming deuteron, that 
of the intermediate two-nucleon system. For our numerical evaluations we use the masses^ 
me = 2mAT, md = 2m^ — B. As we calculate in the cm frame of the yd system, = 0 in 

the s-channel diagram, whereas Pc = —ki — kf in the n-channel, i.e. 2 ^ = ^ (1 + cos 9). 
— Ec denotes the excitation energy of the intermediate state C. 

The interaction Hamiltonian in Eq. (El is 


//“=- / /(O'^K)rfy 


(2.5) 


An explicit expression for the photon field, expanded into multipoles, has been derived 
in d in analogy to Ref. ET! and is given in Eq. dni). It consists of three parts and can 

^Isospin-breaking effects are neglected, i.e. we set nip = nin = mjsf. 
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Figure 6: Diagrams with two-nucleon intermediate states in the s-channel (left) and in the 
M-channel (right). The labels denote the energies of the various particles and states. 


schematically be written as 

1 = V0 + (2.6) 

with denoting the magnetic part of A, V(p + A^'^'> the electric part. The dehnitions of the 
function 0(^), as well as of A^^'> and A^^\ are given in Appendix^ Now we systematically 
replace the photon held in the interaction Hamiltonian (Q by the three terms contained 
in Eq. (ESI). Therefore, when we only distinguish between the various possibilities for A, we 
hnd nine different combinations in Eq. (E3D. The largest contributions are those where we 
substitute A(,^) —^ V0(.^) at both vertices, which is the only part of the photon held that 
contributes for a; = 0, cf. Appendix El Further terms, where this replacement is made only 
once, are discussed in Section 12.2.,‘fl Only a few combinations of interactions without the 
gradient part of A give visible contributions^. These are also taken care of in Section [2. 2. dl 
In this section, we calculate Eq. ()2.4|1 with FT”* ~ f -^(O ' V0(^) simultaneously 
at both vertices, i.e. we restrict ourselves to the terms arising from minimal coupling. In 
order to simplify the calculation on the one hand, and to ensure gauge invariance and the 
correct Thomson limit on the other, we integrate by parts and use current conservation: 


/(O-v0(Od3e = 

j V-j(O0(Oci'e 

(2.7) 

V ./(?) = . 

= [»,«>] 

(2.8) 


The fact that one only needs to know the charge density p in order to calculate the amplitude 
in the long-wavelength limit is referred to as “Siegert’s theorem” inn. For p{^) one can hnd 

^Contributions to the elastic cross sections of the order of < 2% are neglected throughout this work. 
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in [3H] the general decomposition 

p(0 = + (2.9) 

with the charge density of the two nucleons and Xp, Xn the position of proton and neu¬ 
tron, respectively. Xp, Xn) is the charge density arising from meson-exchange currents, 

e.g. from those given in the lower line of Fig. [2l The dominant term in Eq. (j2.9j) is 

= e(5(f-fp), (2.10) 

j=n,p 

which is the only non-vanishing contribution to p {^) in the static limit (“Siegert’s hypothe¬ 
sis” |H5). Note that the ^-functions in Eq. fiTTOl) indicate that the two nucleons are treated 
as pointlike particles, i.e. unlike the authors of e.g. Ref. na we do not introduce any 
nucleon form factors. We also performed calculations including Xp,x„). From these 

investigations, which are reported in Ref. [2H]; we conclude that is well negligible in the 
energy range considered. Therefore we are only concerned with ) in this work. 

Albeit it is not obvious, the use of current conservation in Eq. (UHl) causes that meson- 
exchange currents (cf. Fig. [7j) are also implicitly included in the calculation, as 

V ■ = -t [W" ri ■ T2, - i [H, p^^] . (2.11) 

is the np-potential from one-pion exchange which is part of the Hamiltonian H. fi 
is the isospin operator of the Rh nucleon. These contributions, which go beyond coupling 
the photon held to the single-nucleon current, like in Fig.Q were already indicated in Fig. [21 


2 1 
(a) 


1 2 
(b) 


Figure 7: One-pion-exchange currents contributing to our calculation: the “Kroll-Ruderman 
current” (a) and the “pion-pole current” (b). 


Substituting V ■ J by —i 


in Eg. (12.7|) . the integral over the dummy variable 


^ can easily be performed to yield 


ipmt ^ e 0(xp)], 


( 2 . 12 ) 
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where H is the full Hamiltonian of the np-system 


H 


Pp I Pn 
2mp 2m„ 




(2,13) 


with the np-potential V. In order to evaluate the commutator we switch to cm 

variables, i.e. 


Pp Pn fj ^ ^ ^ Xp Xn /'o 1 /|^ 

P = -^—, P = Pp + Pn, r = Xp-Xn, -• ( 2 - 14 ) 

Our analysis shows that recoil corrections, which arise from the cm motion of the deuteron, 
are well negligible. This observation agrees with Ref. uni, where such corrections have been 
evaluated as well. The net effect is that we may set R = 0, i.e. we neglect the cm velocity of 
the two nucleons'^. As a consequence, we hnd Xp = r{2 and the Hamiltonian ()2.1dj) simplihes 
to the “internal” Hamiltonian 

~*2 

= — + V. (2.15) 

rriN 


Therefore, we can rewrite Eq. (ITT^ as 

H”* = -i [H^P, e 0(r/2)]. (2.16) 

Inserting the commutator (innD into Eq. dZl and dehning cf)i = e0j(r/2), 0/ = e0/(r/2) 
in analogy to Ref. CB] we get 




E 

c 


{df\ \C){C 




d^) 


uj + ^-B-Ec 

2md ^ 


{d,\ |C)(C| |(i.) 


-^ + £:-b-Ec-£s^ 


(2.17) 


In order to keep track of the various combinations of interaction Hamiltonians we have 
labelled the double-^ transition matrix ’00’; the photon states have been skipped for brevity. 
Now the commutators are expanded and, as | dij), \ C) are eigenstates of we can 
act with on these states. We end up with four amplitudes, which have already been 


^We note that there is one exception to this rule: We do include the kinetic energy of the two-nucleon 
system in energy denominators, cf. Eq. 
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derived (in the lab system) in Ref. |lt)j : 




V\^{df\4>f\C){C\cl>i\d,) 


UJ + 


U! + 


? 


2md 

^2 
Pc 

2mc 2md 


c 


UJ 


uj+^-B-Ec 

2md ^ 


E 


{df \^i \ C){C \^f \ di) 


c —UJ 


Pc 


2mc 2md 


B-E, 


c 


Pc 

2md nid 


I d f I 0/ I di ) 


■idj 




, 0 / + 


di) 


(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


is the only one of these amplitudes which contributes in the static limit. It is respon¬ 
sible for the correct low-energy behaviour of the calculation, as will be discussed in detail 
in Section r2. 2. 41 

We defer the evaluation of these amplitudes to Appendix |Bl and turn now to those 
contributions, where the substitution A(^) —>• V 0 (^) (cf. the beginning of this section) is 
made at most once. 


2.2.3 Subleading Terms 


So far we only considered contributions arising from minimal coupling of the photon held 
to the two-nucleon system at both vertices. In the following, we describe how to calculate 
the amplitudes given in Eq. (El, when the replacement 


Hint 


) ■ A(e) d^^ 


/(O-v 0 (o 


( 2 . 22 ) 


is made only once. The term ’subleading’ refers to the fact that the resulting amplitude is 
numerically less important than that of Section 12 .2.21 - its contribution to the differential 
cross section is suppressed with respect to the dominant terms from Section l2.2.2l bv at least 
one order of magnitude for all energies and angles considered, see Fig. El The amplitude 
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is denoted by ■M'fi and follows immediately from Eq. (Q; 




{d,\S /(f) ■ I C) (C I / J(f). A (^) (Pa I d, > 


c 


+ Ik; --B - 


{d, I / /(f) ■ i(f) I c)(c I / j(f) ■ v*(f) d^a 14 > 


{df\ 

oo + ^-B- 

2md 

s ■m-Vd>,(0'H\c){c\ 

Ec 

/J(0-1(0 d^a 

d^) 

{df 

1 / M ) ■ M ) ‘Pi 

Pc 

d 2mc 

C){C 

-B-Ec 

\}j(a)-v,i,f(i)d^( 

\d^) 


-00 + 


_ nc 
2md 2mc 


B-E, 


(2.23) 


c 


Now we perform the same steps as described in Eqs. (I2.12II2.21I1 . i.e. we first replace J J(0 ' 
V0(.f) by i e0(r/2)], then act with H^'p on | d) and | (7), respectively, and hnally 
add and subtract terms in order to perform some cancellations against the denominator. 
We hnd, again neglecting recoil terms and the deuteron velocity, 






j(a)-A(i)d^a\di) 


c 

- ^o; + 
'df\ 
+ ( + 


\ {df I 0; I C){C IJJia-AiOd^^ I d.) 


oo + ^-B-Ec 

2md 


f ■m-A(()d^(\C){C\h\di) 

^^\{d,\ //(?)■ A(t)d«!; I C)(C \4.,\d,) 


2md 




+ (<i,|*|C>{C| / 14) 


^ Pc\{df\<t>.\c){c\fj{a)-A{a)d‘a\di) 

00 — — -h 


2md 2mc 


-00 + 


jE_ _ 

2md 2mc 


B-E, 


c 


{df\ / m-A{i)d^i\c){c\^f\d,) 


00 


oo^ ^ Pc Yd^\j J{i)-A{i)d^i\C){C\ct 0 f\d,) 


2md 2mc J 


oc 


(2.24) 


Whenever the energy denominator has been cancelled, the sum over C may be collapsed. 
As 0 J{^) ■ A(,f) d^^ ~ ) ■ ^('^) 4>i these four terms cancel exactly, and only 

the terms including an energy denominator remain. 
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The relevant parts of the photon held ) are the non-gradient terms in Eq. (12.dj) . The 
cnrrent J(^) includes one-body and two-body currents. The latter are the pion-exchange 
currents of Fig. [71 however we found that only the Kroll-Ruderman current, Fig. (TJ^a), gives 
visible contributions. The single-nucleon current consists of two parts, which we call 
and with 




J=n,p 

j=n,p 

cf. e.g. j3H]. So we may schematically write 


J(f) = J»(f) + J»(f) + /KR(f), 


(2.25) 

(2.26) 


(2.27) 


fij is the magnetic moment, aj the spin operator and p) the momentum of the jth nucleon. 
We observe sizeable contributions only from the magnetic moment interaction encoded in 
the spin current Eq. (I2.25[) . whereas turned out negligibly small. Our notation for 
the amplitudes is when we replace J{^) by ) in Eq. (j2.24jl and when we 

use the Kroll-Ruderman current instead. As this part of our work is rather technical, we 
shift it to Appendix El There, we also discuss those amplitudes which do not contain the 
gradient part of the photon held but nevertheless give sizeable contributions, i.e. Af j/" and 

Two-body currents with explicit A(1232) degrees of freedom, as displayed in Fig. jSl are 
suppressed by one order in e with respect to the Kroll-Ruderman current. Fig. |7Ka), due 
to the 'jNA vertex being part of [H]. This agrees with the hndings of na. where 
such contributions to elastic deuteron Compton scattering below 100 MeV were claimed to 
be of the order of 2 %. A similar size is reported in |2H] for the process np —> dy, where 
the contributions from the A(1232) current turn out to be considerably smaller than those 
from pionic exchange currents. Therefore, and due to the excellent agreement of the total 
deuteron-photodisintegration cross section, extracted from our elastic Compton amplitude, 
with data, cf. Appendix O we so far refrain from including these terms into our calculation. 
Nevertheless, it would be an interesting future task to perform a detailed investigation of 
their size. 

Now we have prepared all ingredients of our deuteron Compton calculation. In the next 
section we demonstrate that it fulhlls the well-known low-energy theorem, Eq. (Q, i.e. 
we obtain the correct low-energy limit within our approach. In Section |21 we present our 
results at non-zero energies and compare them to those from Ref. P and to data. 


2.2.4 Low-Energy (Thomson) Limit 

This work constitutes an alternative hybrid approach to the calculations of Refs. 

which are not applicable for photon energies below 50 MeV. In this section, we prove that 
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Figure 8: Exemplary one-pion-exchange currents with explicit A(1232) degrees of freedom. 


we have indeed removed the limitations of these works at low energies, i.e. we reach the 

correct limit of vanishing photon energy, Eq. El . 

The only non-vanishing amplitudes in the static limit - except for the proton seagull, 

Fig.ISfa) - are the explicit pion-exchange diagrams. Fig. [HI and the double-commutator term, 

~*2 

Eq. ()2.21|1 . This double-commutator involves the internal Hamiltonian -|- V, cf. 

Eq. ()2.15j) . and therefore can be separated into a kinetic and a potential part. Arenhovel 
showed analytically that in the static limit, the potential energy part, using the one-pion- 
exchange potential, cancels exactly the contributions from explicit pion exchange. Fig. El 
[32]. Therefore, the kinetic energy part has to cancel half of the proton seagull. This can 
easily be shown to be true, cf. Ref. uni or j2Hl. Note that the Thomson amplitude (El is 
independent of the deuteron wave function and the np-potential chosen. 

Our numerical evaluation agrees well with the Thomson limit (El, as demonstrated in 
a comparison (Fig. (HI) between the proton Compton cross section, = | of this cross 
section and the deuteron Compton cross section at zero photon energy (we remind the reader 
that ^ cx |Af jjp). The latter two curves are nearly indistinguishable. In the right panel of 
Fig. El we see that the relative error (^)^/ {\%)p — 1 is constant in 6 and less than 2%. 
Therefore it can be accounted for by a constant factor. The main part of this discrepancy is 
due to numerical uncertainties in the normalization of the wave function within our code. 

In this section, we showed that our Green’s-function hybrid approach to deuteron Comp¬ 
ton scattering fulhlls the low-energy theorem and therefore guarantees gauge invariance of 
the calculation. In the next section, we present our results for non-zero photon energies. 
From the good description of all data available we conclude that we have achieved a consis¬ 
tent description of yd scattering for photon energies ranging from 0 MeV up to a; ~ 100 MeV. 


3 Predictions and Discussion 

3.1 Comparison to Previous Work 

In Fig. [ini we show our predictions for the elastic deuteron Compton cross sections, com¬ 
pared to those from Ref. [T] and to the data from Illinois [7j, Lund [S] and SAL [2]. Our 
results are parameter-free, as we use the values obtained from proton Compton scattering 
in [Sj for the two a priori unknown parameters gnj and gus, cf. Fig. at), Obviously we 
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Figure 9: Left panel: Comparison of the proton (black, solid), deuteron (grey, solid) and 
1/4 of the proton (black, dashed) Compton cross section in the static limit. The function 
plotted in the right panel is - 1 - 


have reached our final goal: a (chirally) consistent calculation of elastic deuteron Comp¬ 
ton scattering which describes all existing data reasonably well and also satishes the low- 
energy theorem exactly, cf. Section 12.2.41 If not stated differently, all curves throughout 
this work have been derived using the NNLO chiral wave function from Ref. ra with cutoff 
A = 650 MeV. However, as we demonstrate in Section EH we achieve very similar results 
with other state-of-the-art wave functions. The values we use for the various input parame¬ 
ters are given in Tabled The numbers for the short-distance couplings 5 ^ 117 , are taken 
from the Baldin-constrained fit jS] of and (3m to the proton Compton data, cf. Eq. (Q. 
We use the resulting proton polarizabilities for the neutron analogues as well, as there are 
no isovector contributions up to third order in the SSE scheme. 

There are still minor deviations from the experiments, e.g. our calculation lies slightly 
above the three 49 MeV data from [7] which have been measured at angles below 120°. 
However, this is a feature that our calculation has in common with other approaches which 
also reach the correct static limit, e.g. naiiniEni. At higher energies, the two calculations 
of Ref. PP and of this work approach each other. This is another important cross-check as 
the power counting of the calculations of Ref. PP, as well as of Refs. imiiHi, was designed 
for photon energies of cu ~ 100 MeV. Consequently, both curves in Fig. EUl describe the 
94.2 MeV data from jU] equally well - in fact they nearly lie on top of each other. 

Apart from the total result, we discuss the strength of several contributions separately. 
There are however certain amplitudes which are closely related to each other, e.g. the ki¬ 
netic energy part of the double commutator, Eq. cancels half of the proton seagull in 

the static limit, cf. Section F2. 2. 41 The sum of the potential energy part of the commutator 
and the nine two-body contributions from Fig. El is zero in the limit of vanishing photon 
energy. It stays small in the whole energy range considered in this work, as already ob¬ 
served in Refs. ESI uni. Therefore we do not separate these contributions from each other. 
Nevertheless, there are a few issues worth investigating in more detail: 
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Figure 10: Comparison of our predictions for deuteron Compton scattering (black) with 
those from Ref. PP (grey). In both calculations the numbers from the 2-parameter SSE 
£t to proton data of Ref. jS] are used for the isoscalar polarizabilities. The data are from 
[Zj (circle), [Hj (star) and [Oj (diamond). 


1) The prominent role of the amplitudes cf. Eqs. fj2.18[ l7!T?Hl . which include an 

El-interaction at both vertices. 

2) The importance of the amplitudes and Af , with a denoting the coupling to 

the spin current. 

3) The strength of the amplitudes with the explicit Kroll-Ruderman current at one vertex, 

In the upper two panels of Fig. we investigate the two extreme energies of Fig. these 
contributions are successively added to the remaining terms: the single-nucleon amplitudes 
from Figs. El and EJ the two-nucleon diagrams from Fig. El and the double-commutator 
amplitude cf. Eq. (I2.21j) . The amplitude Af (Eq. (I2.2()|l i is a small correction and 

has been added to the leading amplitudes Af 

Obviously, the amplitudes Al^f are the dominant ones in Fig. El This observation 
holds for both energies considered. However, the amplitudes Af give also important con¬ 
tributions. The same pattern occurs in the calculation of total deuteron-photodisintegration 
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Parameter 

Value 

Comment 


139.6 MeV 

charged pion mass 

niN 

938.9 MeV 

isoscalar nucleon mass 

P 

0.075 

pion-nucleon coupling constant 

fn 

92.4 MeV 

pion-decay constant 

9A 

1.267 

axial coupling constant 

a 

1/137 

QED hne-structure constant 

hp 

2.795 

magnetic moment of the proton 


-1.913 

magnetic moment of the neutron 

Aq 

271.1 MeV 

AA mass splitting 

9ttNA 

1.125 

ttAA coupling constant 

hi 

4.67 

yAA coupling constant 

9in 

18.82 

short-distance coupling constant 

5'118 

-6.05 

short-distance coupling constant 

rrid 

1875.58 MeV 

deuteron mass 

B 

2.2246 MeV 

deuteron binding energy 


Table 1; xEFT parameters determined independently of this work. Magnetic moments are 
given in nnclear magnetons. 

cross sections, cf. Appendix O The contribntions from Af are nearly negligible. The 
small size of these terms is dne to the fact that the two amplitndes which arise from conpling 
the two non-gradient parts of A to the spin current, cf. Eqs. (EH 1123), largely cancel each 
other. The diagrams with one photon explicitly coupling to the Kroll-Ruderman current 
are tiny for low energies, but give a sizeable correction at 94.2 MeV. This contribution is 
stronger in our calculation than it appears in which one may partly attribute to the 
fact that we do not neglect the photon energy in the denominator of the pion propagator of 
the Kroll-Ruderman current, in contradistinction to CB], see Eqs. 

We also give an estimate of the strength of contributions from photons with multipolarity 
L = 2, cf. Eq. ()A.1|1 . In Ref. these next-to-leading terms in the multipole expansion of 
the photon held are claimed to be small and therefore have been neglected. However, we 
slightly disagree from this statement, as can be seen in the lower two panels of Fig. where 
we compare our full results to curves which only include the L = L' = 1-approximation of the 
dominant amplitudes and For low energies, these corrections are certainly 

negligible, but they are of the order of 10% in the high-energy regime in the forward and 
backward direction. 

Comparing to Ref. CB!, we see the main difference to this work in our systematic yEFT 
description of the single-nucleon structure. In CB!, the structure of the nucleon is included 
only via the static polarizabilities and Pm, he. via the leading terms of a Taylor ex¬ 
pansion of the single-nucleon Compton amplitudes. In our work, these amplitudes have 
been calculated up to third order in the Small Scale Expansion, as explained in detail in 
Refs. [21 C] , and are included with their full energy dependence. Another advantage of our 
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Figure 11: Comparison of separate contributions to our deuteron Compton-scattering re¬ 
sults. In the upper panels we compare the full result (solid) to curves with several ampli¬ 
tudes subtracted; the subtracted amplitudes are: Al™ (dotdashed), (dotted), 

_/VfKR _|_ (longdashed), A4f^ + (shortdashed). In the 

lower panels, we compare our full result (black) to a curve where the amplitudes Mfl and 
have been replaced by their L = L' = 1-approximations (grey). The data are from 
[7j (circle) and [H] (diamond). 


approach is the treatment of the pion propagator in the pion-exchange diagrams of Fig. 0 
We calculate these diagrams using the full pion propagator, whereas the authors of M 
always assume that the photon energy is small compared to the energy of the virtual pion 
and therefore may be neglected. This, however, is no longer a good approximation as the 
photon energy approaches the pion mass. A similar difference occurs in the Kroll-Ruderman 
amplitudes, as discussed above. Finally, we do not agree with the statement of M that 
L = 2-contributions are negligible for all amplitudes and energies considered, see Fisr. 1111 
We showed in the last two sections that our calculation provides a consistent description 
of elastic deuteron Compton scattering below 100 MeV. It will be interesting to see whether 
data from a forthcoming experiment at MAXlab at a; ~ 120 MeV are closer to the Green’s- 
function hybrid approach or Weinberg’s hybrid approach. There, however, corrections due 
to the kinematically correct position of the pion-production threshold in analogy to Ref. [3], 
which are not yet included in our calculation, should not be neglected. Therefore, the 
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120 MeV curve in the direct comparison of our results at various energies in Fig. is only 
a qualitative statement about the behaviour of the differential cross section for uj —> m^. 
Hildebrandt et ah pp estimate this effect to be negligible below 100 MeV. We also show 
our prediction at cuiab = 30 MeV, which might be considered as a guideline to forthcoming 
experiments in this energy region HH. It is comparable in magnitude to the 49 MeV curve. 
Its shape, however, is less asymmetric between the forward and the backward direction, due 
to the exact forward-backward symmetry of the static cross section, cf. Fig. 1^1 



^ab [deg] 

Figure 12: Comparison of our results for differential Compton cross sections at various 
energies: 30 MeV (grey), 49 MeV (shortdashed), 68 MeV (dotted), 94.2 MeV (dotdashed), 
120 MeV (longdashed). 


In the next two sections, we investigate the sensitivity of our results to the deuteron 
wave function and the W-potential used. 

3.2 Dependence on the Deuteron Wave Function 

As demonstrated in Section r2. 2. 41 our calculation fulhlls the low-energy theorem, Eq. (I2.2|l . 
which in turn is independent of the deuteron wave function. Therefore and because of the 
nearly energy-independent offset between the cross sections calculated with the chiral ra 
and the AV18-wave function observed in Refs. IIII2H1, it is not surprising that the 
wave-function dependence of our present calculation is also at non-zero energies considerably 
reduced with respect to Refs. EIIHI. In fact, the remaining dependence is of the order of 
1% and therefore nearly invisible, cf. Fig. where we compare our cross sections with 
the two wave functions that turned out as the extreme ones in Ref. : the AV18 Q 
and the NNLO yPT [2^] wave function (the same observation holds for other state-of-the- 
art deuteron wave functions). This is another important success of our present approach 
to deuteron Compton scattering, as it demonstrates that our calculation is not sensitive 
to details of high-energy physics, i.e. short-distance contributions of the wave function. 
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whereas the 10%-effect observed in Refs. manifests a much stronger dependence on 

short-distance dynamics than would be expected from a low-energy Effective Field Theory. 
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Figure 13: Comparison of our deuteron Compton cross-section results for 68 and 94.2 MeV, 
using two different wave functions: NNLO yPT (grey) [HEl and AV18 (dashed) j22j. In the 
lower two panels we show / (f )avi8 “ I’ 


3.3 Dependence on the Potential 

In this section, we investigate briefly the sensitivity to the np-potential used for the rescat¬ 
tering part. Usually, we use the AVI 8-potential ra which provides an excellent theoretical 
description of the Nijmegen partial-wave analysis. We compare here our results achieved 
with this modern ’high-precision’ potential to those of the leading-order chiral potential, 
which includes only the one-pion exchange and a simple parameterization of short-distance 
effects via two point-like, momentum-independent contact operators. This potential is given 
e.g. in Ref. EOl, using a Gaussian regulator in order to render the pion-exchange potential 
hnite at the origin. 

At leading order, there are two free parameters, Cif° and Cq, with d denoting the 
(deuteron) ^Si-^Di channel. has been hxed in via the ^S'o scattering length, oq ~ 

—23.75 fm, Cq to the deuteron binding energy, parameterizing any short-distance physics 
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in the spin-singlet and -triplet channel, respectively. Cq^° and Cq, as given in Ref. [101 for 
the cutoff-value A = 600 MeV, are reported in Table |21 


^0 


-0.422 fm 

0.795 fm 


Table 2: Parameters of the LO chiral potential as determined in jlU] for A = 600 MeV. 

Even with this rather crude approximation of the neutron-proton interaction, we obtain 
results close to those of the AV18-potential, cf. Fig. El Obviously, the one-pion-exchange 
potential, adequately regulated for r —> 0 , together with a reasonable parameterization of 
the hard core gives an approximation of the potential which is well sufficient for the process 
under consideration. We conclude that we are mainly sensitive to the long-range part of 
the potential. Nevertheless, there are small deviations (of the order of < 4%) visible in 
Fig. HH which justify the application of a more sophisticated potential. Not surprisingly, 
these deviations, which arise due to the poor description of high-energy (short-distance) 
dynamics in the LO chiral potential, increase with increasing photon energy and scattering 
angle. 

Having proven that our calculation is rather insensitive to the choice of the deuteron wave 
function and the np-potential, we turn now to hts of the isoscalar nucleon polarizabilities (or 
equivalently the short-distance contributions 5^117 and pus) to all existing elastic ■yd data. 


4 Fits of the Isoscalar Polarizabilities 

We saw in Section ED that our results for the elastic deuteron Compton cross sections 
give a good description of all existing data. Furthermore, the cross-check described in 
Appendix O i-e. extracting the total deuteron-photodisintegration cross section from our 
Compton amplitude, together with the exact reproduction of the low-energy theorem, cf. 
Section mu gives a strong hint that the numerically most important amplitudes have 
been calculated correctly. Therefore, we now use our deuteron Compton cross sections 
to fit the static isoscalar nucleon polarizabilities and to elastic yd experiments. 
This corresponds to htting the coupling constants 5^117 and gns of the two short-distance 
operators. Fig. Hf). to elastic yd rather than yp data. We can use all data for the hts, 
whereas the authors of Ref. P] had to restrict themselves to the experiments performed 
around 68 and 94.2 MeV and those of Ref. US] excluded the two 94.2 MeV data in the 
backward direction in certain hts. 

We do a least-y;^ ht, using the chiral NNLO wave function with A = 650 MeV EHI- Our 
results for the isoscalar polarizabilities from the global ht to all data read 

= (11.5 ± 1.4 (stat)) ■ lO”"^ fm^, 

Pm = (3-4 ± 1.6 (stat)) ■ 10“^ fm^. (4.1) 
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Figure 14: Upper panels: Comparison of our result at 68 MeV and 94.2 MeV using two 
different np-potentials: the AV18-potential 1221 (solid) and the LO chiral potential UDI 
(dashed). For both curves the chiral wave function has been used. Lower panels: 
Corresponding error plots, calculated in analogy to Fig. 121 


We only give the statistical error as we neglect further uncertainties, e.g. the error induced 
by the dependence on the deuteron wave function. This error may well be set to zero, 
due to the tiny wave-function dependence observed in Fig. whereas the wave function 
introduces a sizeable uncertainty in the calculations of Refs. PIIHI. Theoretical errors from 
higher orders are also neglected, as we could only give a rough estimate of their magnitude®. 
Nevertheless, we are aware that they may be comparable in size with our statistical error. 
The corresponding per degree of freedom is 

= 0.98 (4.2) 

d.o.f. 

with 27 degrees of freedom (4 data points from [7] at 49 MeV, 9 from [S] at 55 MeV, 2 from 
[2j, 9 from [Hj around 68 MeV and 5 from [2j at 94.2 MeV, along with two £t parameters). 
In Fig. Uni we give a contour plot of the achieved together with the 70%-conhdence 

®In pp, such corrections have been estimated to be of the order of 1 • 10“^ fm"^ from naive dimensional 
analysis. 
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Figure 15: Contour plot of the achieved in our global 2-parameter fit for and /3f^ 
(left panel) and the 70%-confidence ellipse (right panel). The shadings in the left panel are 
separated by steps of = 0.67. 


The plots corresponding to Eq. m, together with the (statistical) error bands, are 
shown in Fig. The predictions within our Green’s-function hybrid approach, using 
the results for and from Eq. (11.211 for the proton and the neutron polarizabilities, 
describe the data already well, see Fig. (TUI It is therefore no surprise that also the fitted 
curves are in good agreement with experiment. We compare our fit results to “fit IV” from 
Ref. [IS], which is the 0{q^) HB^PT fit to all data with central values = 11.5 ■ 10“^ fm^, 

= 0.3 ■ 10“^ fm^. As explained in detail in PP, the only sizeable deviations are observed 
at 94.2 MeV in the backward direction, due to the A-resonance diagram. Fig. HKa), which 
is not included in the calculation of Ref. ^Hj- In Fig- Eland Table |S1 we also compare to 
the 2-parameter ht from Ref. P], which was performed with the yd —> yd kernel according 
to third-order SSE, using the chiral wave function jSSj^- Here we observe a constant offset 
in the differential cross section at 68 MeV, whereas at 94.2 MeV the two curves are quite 
close to each other, similarly to Fig. El 

The value of our global fit for a% is slightly smaller, the one for slightly larger than 
the fit results of Ref. PP, given in Eq. (Q and in Table 0] respectively. Nevertheless, both 
extractions agree well with each other within their error bars, and there is also very good 
agreement of Eq. (03) with the values given in |I|, see Eqs. 03 and 031 ). Furthermore, 
we hnd that the numbers given in Eq. 03 add up nearly exactly to the isoscalar Baldin 
sum rule (see PP for the proton and neutron sums used as input). 


+ ( 3 ^ 


M 


= (14.5 ± 0.6) ■ 10“^ fm 


world av. 


(4.3) 


which has been a serious problem in former extractions O Ej. Therefore, in order to 
reduce the statistical error, we repeat our global ht, using the central sum-rule value as an 

®We note that these curves correspond to the fits performed with the “effective” data set, cf. Table 3 
and Fig. 11 in p. 
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Figure 16: Results from a global fit of a|; and to all existing elastic 'yd data (solid). 
The grey bands are derived from our statistical errors. The dotted line represents “ht IV”, 
one of the (9(g'^)-HBxPT hts from Ref. [TH], with central values a% = 11.5 ■ 10“^ fm^, = 
0.3 ■ 10“^ fm^. For the 0{q'^) calculation the NLO chiral wave function of Ref. [H] has been 
used, whereas our results were derived with the NNLO-version of this wave function |36j . 


additional £t constraint. The results are 


cy. I 


A 


M 


Baldin 


Baldin 


(11.3 ± 0.7 (stat) ± 0.6 (Baldin)) ■ 10 fm^, 
(3.2 =F 0.7 (stat) ± 0.6 (Baldin)) ■ lO”'^ fm^ 


(4.4) 


with /d.o.f. = 0.95. Of course, the central values of Eq. (I4.4|l are very similar to the ones of 
Eq. (14.111 . due to the nearly perfect agreement of the 2-parameter-ht result with the sum-rule 
value. However, the statistical error of the Baldin-constrained fit is reduced by about 50%. 
The plots arising from the global, Baldin-constrained ht, together with the corresponding 
error bars, are shown in Fig. The central curves are nearly indistinguishable from 
the ones of Fig. [TBl In order to simplify comparison, we sum up our fit results, together 
with those from Refs. in Table El Obviously, all three extractions of the electric 

polarizability agree with each other, whereas determining I3m without the explicit inclusion 
of A resonance degrees of freedom yields a negative central value in contradiction 
to |31 E] and the present work, which show also good agreement in this quantity. 
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Figure 17: Results from a global fit of a% and to all existing elastic 'jd data (solid), 
obtained with the NNLO chiral wave function The grey bands are derived from our 
statistical errors. We compare here our results to the 2-parameter-£t results from P] , using 
the same chiral wave function (dashed). 



a% [10“"^ fm^] 

Hi, [10- fm’] 

III 

12.8 ± 1.4 (stat) ± 1.1 (syst) 

2.1 ± 1.7 (stat) ± 0.1 (syst) 

nBi 

13.0 ± 1.9 (stat)+i ;5 (syst) 

—1.8 ± 1.9 (stat)lo !9 (syst) 

this work 

11.5 ± 1.4 (stat) 

3.4 ± 1.6 (stat) 


Table 3: Comparison of the fit results for the isoscalar nucleon polarizabilities a% and 
achieved within 0{e^) SSE p[] (Weinberg hybrid approach), 0{p^) HByPT ^H] (Weinberg 
hybrid approach) and the Green’s-function hybrid approach presented in this work, respec¬ 
tively. ’syst’ denotes the systematic uncertainty from the wave-function dependence EIIHI 
and the arbitrariness as to which data are included in the ht [TH] . 


Combining Eqs. (M or (D, respectively, with the (Baldin-constrained) fit results of 
Eq. (I1.2jl . which we extracted from proton data using the SSE framework, we calculate the 
values for the neutron polarizabilities as 

Og = (12.0 ± 2.0 (stat) ± 0.4 (Baldin)) ■ 10“^ fm^, 

= (4.0 ± 2.1 (stat) ± 0.4 (Baldin)) ■ 10“^ fm^ (4.5) 


for the 2-parameter £t and 


a, 


Baldin 


= (11.6 ± 1.5 (stat) ± 0.6 (Baldin)) ■ 10 ^ fm"^. 


p: 


M 


Baldin 


= (3.6 =F 1.5 (stat) ± 0.6 (Baldin)) ■ 10 fm' 


(4.6) 


for the fit of the isoscalar polarizabilities including the Baldin constraint. We consider these 
values to be as reliable as those from the quasi-free Compton experiment of Ref. dl, as the 

















Figure 18: Results from a global fit of to all existing elastic 'yd data, using the chiral 
wave function j2nj. is hxed via the Baldin sum rule, Eq. ()4.8j) . The grey bands are 
derived from our statistical errors. 


isoscalar polarizabilities, from which they are derived, have been determined by htting our 
deuteron Compton calculation, which fulhlls the low-energy theorem, to all existing elastic 
deuteron Compton-scattering data. That means there is no restriction on either the energy, 
as in Ref. P], or on the angle, like in Ref. ^H]- From these values we deduce that the 
magnetic response of the neutron is comparable to that of the proton and that nucleons are 
paramagnetic. We also conclude that the isovector polarizabilities are considerably smaller 
than the isoscalar ones. In other words, our analysis shows that within the precision of the 
currently existing data, elastic Compton scattering from the proton and the deuteron is in 
agreement with 

~ Og, ~ /^M- (4-7) 

These findings agree with those of Refs. min. 

5 Conclusion 

In this work we examined elastic Compton scattering from the deuteron. Diagrams without 
an intermediate two-nucleon state have been calculated up to next-to-leading order in the 
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Small Scale Expansion, an Effective Field Theory with nucleons, pions and the A(1232) 
resonance as explicit degrees of freedom. Those diagrams including the propagation of the 
two-nucleon system between the two photon vertices have been calculated using Green’s- 
function methods. Therefore, and because of the fact that we use deuteron wave functions 
that have been derived from state-of-the-art iViV-potentials, we refer to this approach as 
the “Green’s-function hybrid approach”. For the photon coupling we make use of Siegert’s 
theorem [31], which is well-known to guarantee the exact static limit, see e.g. [iniE21E2|. 

In Section im we show that we have achieved a consistent description of elastic deuteron 
Gompton scattering which is valid from 0 MeV up to cn ~ 100 MeV. The advantage of our 
calculation with respect to Refs. PIITIIIHI is that we obtain a good description of the data 
published below 60 MeV. This improvement is of course connected to the correct static limit, 
which had not been obtained in those publications. Therefore, other calculations reaching 
this limit and thus being manifestly gauge invariant are also able to describe the low-energy 
data well, see e.g. [isiiiniisni. However, we achieve good agreement also with high-energy 
data, measured around 94.2 MeV [Hj. In this energy regime the work of Ref. HE! fails and the 
“pion-less” Effective Field Theory used in jEEj is inapplicable. The authors of PI have also 
problems to describe the data in the backward direction without introducing surprisingly 
large isovector polarizabilities. It was already shown in pp that the main difference between 
their approach and ours is in the energy dependence of the resonant Gompton multipoles, 
which is well captured in our calculation due to the inclusion of explicit A-resonance degrees 
of freedom, whereas the authors of Ref. na include the polarizabilities only via the leading 
and subleading terms of a Taylor expansion in the photon energy. 

Having achieved a good description of all elastic deuteron Gompton-scattering data 
enabled us to perform a global £t of the isoscalar polarizabilities to all existing data points, 
published in muM Our 2-parameter-£t results are (see Table ED 

a% = (11.5 ± 1.4 (stat)) ■ 10“'^ fm^. 

Pm — (3-4 ± 1.6 (stat)) ■ 10“^ fm^. (5.1) 


We give only the statistical error as we cannot quantify the theoretical uncertainty from 
higher orders in our Green’s-function hybrid approach and because further uncertainties, 
like the wave-function dependence observed in Refs. iniiHj, are tiny in this approach, as 
shown in Sect ions 13.21 and Id. 31 We note that the numbers of Eq. (lEID are very close to the ht 
results da for the proton polarizabilities, determined in |3| within the SSE framework. This 
leaves little space for large isovector polarizabilities. Furthermore they agree extraordinarily 
well with the isoscalar Baldin-sum-rule value, Eq. dOD- Therefore, in order to reduce the 
statistical error, we repeated our hts including this constraint: 




Baldin 


= (11.3 ± 0.7 (stat) ± 0.6 (Baldin)) ■ 10 fnv. 


P\ 


M 


Baldin 


= (3.2 =F 0.7 (stat) ± 0.6 (Baldin)) ■ 10 ^ fm'’. 


(5.2) 


Gombining Eqs. dEID or respectively, with the proton numbers of Eq. (O, we ob- 
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tained the neutron polarizabilities 


= (12.0 ± 2.0 (stat) ± 0.4 (Baldin)) ■ 10 fni^, 
Pm = (4-0 ± 2.1 (stat) ± 0.4 (Baldin)) ■ 10“^ 


for the 2 -parameter £t and 


al 


p: 


M 


Baldin 


Baldin 


(11.6 ± 1.5 (stat) ± 0.6 (Baldin)) ■ 10 fm^, 
(3.6 =F 1.5 (stat) ± 0.6 (Baldin)) ■ 10“'^ fm^ 


(5.3) 


(5.4) 


for the 1-parameter £t including the Baldin constraint. In both hts the fact that our deuteron 
Compton-scattering calculation is applicable in the whole energy range from 0 MeV to 
100 MeV enables us to include all data into our fit of the isoscalar polarizabilities. From these 
results we deduce that the neutron is paramagnetic and that proton and neutron behave 
rather similar when exposed to external electromagnetic fields. In both points our deuteron 
Compton calculation agrees with Refs. gunnn and with the Chiral Perturbation Theory 
prediction that contributions to the isovector polarizabilities only start beyond leading-one- 
loop order. 

Finally, we strongly advocate enlarging the data base for elastic Compton scattering 
on the deuteron. If further experiments, as planned at TUNL/HI 7 S, at the S-DALINAC 
and at MAXlab, provide additional data below the pion threshold, the increased statistics 
would reduce the statistical error in our fit of the isoscalar polarizabilities. Our results for 
aE and Pm are in excellent agreement with those from the “global” analysis j3] which did, 
however, not include elastic deuteron Compton scattering due to the alleged discrepancy 
between theory and elastic yd experiments. This discrepancy has now been resolved, both 
in P and in the Green’s-function hybrid approach presented in this work. Clearly, reduced 
error bars in coherent Compton scattering from the deuteron would enable us to quantify 
by how much proton and neutron polarizabilities differ. 
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A Multipole Expansion of the Photon Field 

For the photon held ^(0; ^^e multipole expansion derived in in analogy to EH, 

see also Ref. EHl- The result is 


_ 


L =1 M=-L 


■L 


'27r ( 2 L + 1 ) 
T(L + 1 ) 


(A.l) 


X — -*^ejLK)TLM(e) -XLYlm{OJl{ooO • 


u 


€\ denotes the photon polarization vector in the spherical (helicity) basis. The functions 
Wigner D-functions for a = 7 = 0, Flm denote the spherical 

harmonics and 'ipL^uir) = (1 + r -^) jL^ur) with the spherical Bessel functions Jl^z) = 
(- 2 ) and Jiiz) the Bessel functions of the hrst kind. In the static (long-wavelength) 
limit only the gradient term in Eq. (El) survives, as for a; —> 0 , ji(a;r) —> |a;r and 
'0i(a;r) —> This term turns out to be the dominant part of the photon held for all 


energies under consideration. Therefore, like in Ref. [1^], we dehne two scalar functions 

A /27r ( 2 L + 1 ) ^ 

(Pi[n = -2_^ 2_^ 0M,Xi^—\l TtT , IN ^L[uJr)YLM[r), 


L =1 M=-L 
00 U 


00 y L (L -|- 1 ) 




L'=l M'=-L' 


which allow us to write 


eA,:e 




00 y U [U + 1 ) 


v*(f), 


7 «W/( 0 . 


(A.2) 


(A,3) 


We want to decompose the photon held in its electric and magnetic part. Therefore, we 
write Eq. (El) in Eq. (12. bB schematically as A = V(/) - 1 - with 


Y ^ \l L{L + 1 ) 

X [%x ^M,x - 4,;^ (-1)^+^ dM,-x{^) 

00 L 

= - YYY 


2n + 


(A.4) 


X 


L(L + 1 ) 

,7,<5M,A-a4(-l)^+"d4_,(0) 


27r (2L -h 1) ^L+i 


k,X=±l ■^'=1 fo[=-L 

d /■ 1^L+A .L 


(A.5) 


constitutes the magnetic part of the photon held, V 0 -|- is the electric held ra. 
The operators al ^ {<^kx) create (destroy) a photon with momentum k and polarization A. 
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B Calculation of Diagrams with Intermediate 
A^A^-Scattering 


Here we give details about the calculation of the diagrams sketched in Fig. IHl Special em¬ 
phasis is put on the construction of the two-nucleon Green’s function, drawing substantially 

from Ref. uni- 

(Eq. ()2.18j) l will be calculated hrst, in analogy to Ref. [TT)] . Dehning the shortcut 
2 

Eq = uj + — B and neglecting prefactors for the moment we can write this amplitude as 


M 


001 

fi 


OC 


E 

c 


{df I V’l' '^l'm' \ C){C \ 'iPlYlm \ di) 


En — E, 


c 


(B.l) 


where we suppressed the sums over L, M and E, M' for brevity, cf. Eq. (1121). Each wave 
function can be separated into a radial part, denoted by the index ’rad’, and an angular 
part, denoted by a hat. Furthermore, | G) is an eigenstate to with eigenvalue Ec, so 
we can write 

M*f cx {dfraddf I ^ jjnp I Gj-ad G ) ( Grad G | M \ di^^.p^di') . (B. 2 ) 

Crad C ° C 


G is used as a shorthand notation for all angular quantum numbers of the intermediate 
state, i.e. | G) =| Lc Sc Jc Me)- Separating the radial from the angular part of | G) and 
inserting two complete sets of radial states | r) and | r') we get 




E 

Gad C ■ 


X ( r' 


1*00 poo 

0 Jo 

1 


Eq - 


Edr r'‘^dr'{ df 
Grad )( Grad 


c 


Yl' M' I G ) ( d/rad I \ r' ) 
r){r I 'ipL I djrad)(G | Ylm 



(B.3) 


For the deuteron (J = 1) wave function, we use the notation = X]i=o 2 

with uoi^r) {u 2 {r)) denoting the usual radial wave functions u{r) (w{r)) in position space, 
cf. e.g. Ref. jSH]. The indices of the angular wave functions y are /II for orbital angular 
momentum, spin and total angular momentum of the deuteron state, m G {—1, 0,1}. Now 
we can write Eq. dEl as 


M^ff OC EEE {I'llMf I Yl'm' I G)(G I Ylm \lllM,) (B.4) 

Q l=0,2l'=0,2 

C f lot' 1 LOT 

X jj rdrYdYup{Y)'ijjL'{ — ){r'\ _ ^np \ r) Uiir), 

where we have removed the sum over Grad- Integrals without limits are always integrated 
from 0 to inhnity throughout this work. 
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We now have to evaluate the double integral in Eq. (IB.4jl . including the Green’s function 

G^(r, r'; Eq) = {r'\ ^ \r). (B.5) 

However, we need to evaluate the integral for arbitrary functions of r. Therefore we describe 
how to calculate 

= JJ rdrr'dr'ui'{r') Jf{r')G^{r,r';EQ) Ji{r)ui{r). (B.6) 

We do so in two steps and dehne 

= J r'dr'ui^r') Jf{r')G^{r,r';Eo). (B.7) 

Once we have solved this hrst part, it is easy to numerically calculate the remaining integral 


= J rdrui{r) Ji{r) (B.8) 

In order to hnd the function - in the following we use the abbreviation Xci"^) 

brevity - we hrst note that 

(b„ - H"”) Gc(r, r'i £:„) = { r' | r) = (b.9) 

Eq. ()B.9jl dehnes the Green’s function corresponding to Schrodinger’s equation with a central 
potential and the Hamiltonian 




p" 


TUn 


+ VM 


rriM \r dr^ j 


(B.IO) 


The dependence of the potential on the quantum numbers of the interim state | G) is shown 
explicitly. Note that the neutron-proton potential contains a tensor part and therefore not 
only depends on the distance r but on the vector f. The tensor force mixes e.g. the 
deuteron s- and d-states in Schrodinger’s equation. Nevertheless, on the level of the Green’s 
function this matrix equation decouples, cf. Ref. jHHj- The decoupling of Eq. (jEl) guarantees 
that only the diagonal terms of the tensor force contribute. Therefore, the orbital angular 
momentum is well dehned, which allows us to replace L^^Lc{Lc + l) in Eq. 

Eqs. ()B.9fl and (jB.lOfl combine to 


1 Lc {Lc + 1 ) 


niN dr"^ rriNr'^ 


Vc{r) 


rGc{r,r';Eo) = 


6{r' - r) 


which e.g. in the deuteron case reduces to the two differential equations 


1 d^ 
dr"^ 


hrfint (^ 


rGo(r,r';Eo) = 


5{r' - r) 


(B.ll) 


(B.12) 
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(B^13) 



1 £ 

niN dr"^ 


6 

rriN r'^ 



rG 2 {r, r';Eo) 


d{r' - r) 
r 


Kent(’") denotes the central part of the potential, Vten(’") the tensor potential. The indices 
of the Green’s functions in Eqs. f)B.121 IB.13fl reflect the orbital angular momentum state, 
whereas J = 1, S' = 1 is not written down explicitly. Acting with the operator given in 
square brackets in Eq. (IBTTI) on the integral over r' collapses and we hnd 


dr"^ 


+ ruN (Eo - 


Vc(r)) 


Lc {Lc + 1 ) 


rXci^) 


rriN ui>{r) J/(r). 


(B.14) 


This is a second-order differential equation in r with an inhomogeneity, which can be in¬ 
terpreted as a source term. Its solutions are real for Eq < 0 and complex for Eq > 0. 
The latter case corresponds to a; > i?, i.e. the photon carries enough energy to break up 
the deuteron into its two constituents. Obviously, an imaginary part only appears in the 
s-channel diagrams, where the incoming photon is absorbed before the other one is emit¬ 
ted. In Appendix O we will use the imaginary part of the amplitudes to derive the total 
deuteron-photodisintegration cross section via the optical theorem. 

For r —> cx), ui>{r)Jf{r) —> 0 due to Ui>{r) 0, i.e. Eq. (IB.ldfl reduces to a homogeneous 
differential equation. Furthermore, V^{r) —»■ 0 for r —> oo. Therefore, we are for large 
distances left with 


' d^ r. 

— + EoiriN 
dr^ 


Lc {Lc + 1 ) 


rxc(’') = 0, 


(B.15) 


This equation is known to be solved by a linear combination of the spherical Bessel functions 
of the hrst and second kind, jiciQ^) ^^d nL^{Qr) with Q = y/rriN Eq. Note that Q can 
be real or imaginary^, depending on Eq. In our case the boundary condition is that Xci"^) 
must be an outgoing spherical wave for large r, cf. e.g. jl^I- Therefore we may write 

Urn Xcir) h^HiQr), (B.16) 

r^oo ^ 


with {Qr) the spherical Hankel function of the hrst kind, dehned as 

h^liQr) = jLciQr) + iriLciQr). (B.17) 

Once we have numerical solutions for the homogeneous and the inhomogeneous differ¬ 
ential equation, we need to hnd the correct linear combination which satishes the condi¬ 
tion ()B.16jl . In other words we have to determine the coefficient A which fulhlls 

™ oc h^HiQr), (B.18) 

^In order to fix the sign of the imaginary solution, one has to add an infinitesimal imaginary part to —B, 
i.e. B ^ B — ie. 
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where (x| 5 °™(^)) denote the solution to the inhomogeneous (homogeneous) differential 

equation. In the asymptotic limit, must be a linear combination of jLciQ^) ^md 

nLciQ"^) equivalently, of jiciQ^) and h^l^^{Qr). Therefore we can write the general 
solutions in the following way: 


with functions 

o 

With the choice A = 


horn 


XK“(r) = Ca “hr 


■vhom ( 

'c 


jicW’-) + «e”MhcW’-) 

X”(r) = C“(r) iLc(Qr) + y(r) h'-ll(Qr 


(r), which become the constants 

we hnd 

O ' o 


(B.19) 

(B.20) 


,in/hom c i 

tor large r. 
o 


lirn {x”(r) + AxrW}=C” ((“ - t 


hom\ 
C ) 


/‘Ll(Qr) 


(B.21) 


which satishes the condition (IB. 1 fill . Therefore we need to determine the coefficients C^, 

o 

^hom_ 'jpjg pg clone in the region where C^(r), t(^(r) are constant, i.e. their derivatives 

vanish. In this region 


1 in/hom/ \ i /^in/hom , n 

hiXA (r) = lnOA +m 
o o 


jLc(Qr) + t‘"'“”"h'j‘l(Qr) 


(B.22) 


Dehning 


/ in/hom 




c 


n/h( 


in/hom 


we hnd 


£)in/hom _ 

dr 




inJQr) + tf--A/,W(Qr) 


jLc(Qr) + f 


in/hom ^(1) 


C 


Lc 


{Qr) 


(B.23) 


From this equation we can easily determine which we use to solve Eqs. (Ibt^ib:^ 

for and thus to determine A. 

o 

Numerically, this is one of the most involved parts of this work. Fortunately, a nice 
and valuable cross-check to the routine can be performed. For this we consider again the 
double integral to be calculated, Eq. (jB.fif) . This integral is obviously invariant under the 
interchange r <->■ r'. A general feature of Green’s functions is that they are symmetric under 
r r', i.e. G^{r',r; Eq) = G^{r,r'; Eq). Therefore, 


q-ll'C 


r'dr' rdr ui/{r) Jf{r) G^{r, r'] Eq) Ji{r') ui{r' 


(B.24) 


This expression is identical to with i ^ f, I ^ V , i.e. our results must be symmetric 
under i ^ f, I ^ V. This is a non-trivial check, because for J/(r) 7 ^ A(^) completely 
different functions X/^i'^) generated. Our routine agrees well with this symmetry - the 
deviation caused by numerical uncertainties is less than 1 %. 

Now all tools to calculate are prepared. However, as the algebraic manipulations 

are not too complicated, we refer the interested reader to Ref. j2Hj for further evaluation 
and the analytic results. There we also compute and 
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We turn now to the evaluation of the subleading terms of Section \‘2.2.‘Al where we make 
use of the continuity equation (EHD at only one or even none of the two vertices. Therefore, 
we have to specify the current J(^) and the relevant parts of the photon held ). These 
are the non-gradient terms in Eq. (ESD, i.e. (Eq. (lA.dfl i and (Eq. (IA.5jl i. The one- 
body current is considered hrst. It consists of two parts, and J^\ cf. Eqs. fl2.25l 
All possible combinations of and J^'^\ have been calculated in HE]. We also 

evaluated all these amplitudes, however we found that only jA) gives visible contributions to 
the deuteron Compton cross sections. Therefore, we may restrict ourselves to the following 
combinations: {JA'>, denoted by al, and {JA\ A^^), denoted by a2. 

We now calculate J ) ■ A(^’^)(^) We start with the derivation for A^^\ writing 

only the ^-dependent terms for simplicity. 


j=n,p 

= / E [jdujOLYLMii)) (B.25) 


J=n,p 


where one partial integration has been performed. Now we evaluate the integral and after¬ 
wards replace Tp —> |, cf. Eq. (EH, yielding 


j(A(^-')./(1)(^-')^3^0,2[V. X (^j,^{‘^)LYLM{r)')\ (pp Up - (-1) Vn <?n) , (B.26) 


where we used YLMi—r) = (—1)^ By the help of the identity TiLM^r) = , ^ 

v ^ (-f'+l) 

LYLM{f)-i see e.g. jH, this becomes 


JA\^) ■ A(^)(f) <X.2^jL{L + l) Vr jL{‘^)TLLM{r) (pp Up - (-1)^ Un) . 


Now we can use the curl formula jl2l for simplifying 
sion relations for spherical Bessel functions to write 


VrXjL{f)fLLM{f) 


(B.27) 
and the recur- 




■ . /L(L + 1 ) 2 ^ 

2L+1 


. , ujr /L2(L + 1) - 

tuJ3L+i{Y)\/ 2L + 1 


{a. 


P ^p 


(B.28) 


■ 1 ) ^'n ^n) ■ 


We found, like the authors of Ref. ^E] , that the numerical importance of the various contri¬ 
butions rapidly decreases with increasing multipolarity L. Therefore, the term proportional 
to may be neglected. Dehning S = 2>±£ii and t = 22^^ and including all 
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prefactors, we get the result 


J2 Y1 A V2ir (L + 1) ^ i^+i 

fc,A=±l i=l A^=-L ^ 


(B.29) 


X 


(hp ( 1 ) hn) 5'+ (p.p + ( 1 ) fin) t ■ Clk^x^M,\ ;^ ( 1 ) 


The scalar products are replaced using the relation 

fjLMif)-V=[YL<^V]j^, 


(B.30) 


which holds for any vector (rank 1) operator (® denotes the irreducible tensor product). 
An explicit proof of the relation (jB.dnjl is given e.g. in Ref. j2H]. We use it to hnally rewrite 
Eq. (jB.29j) as 


/oo L 

E li E ^V27r(L + l) 

fc,A=±l -^=1 M=-L 


euj 


R+i 


,ujr. 


Itti^ 1 


X { (hp ( 1) hn) [W-l ® ^\lM (hp T ( 1) hn) \Xl-1 ® ^\lM } 


X 




k,\ 


(B.31) 


We turn now to the calculation of J Again we restrict ourselves in 

the derivation to the .^-dependent terms, Ending 

j / E ijL{^oyLM(i)d^( 

j=n,p 

= Vr X [r jL{ — )YLM{f)j {fipdpY {-l)^ findn) , (B.32) 

where we have performed the same steps as in the derivation of Eq. (jESEl). Using the 
relation 


rjL{—)yLM{f) = x/L (L + l)rjL(y) 


ThL-lMif) 


T, 


ll+imY ) 


y(LTlHTLTl) x/L {2L + 1) 
which is derived e.g. in j2H], and the curl formula we find 

TlL-lMi^) TLL+iM{f) 


, (B.33) 


Vr X I rji(—) 


Combining Eqs. 


v/(T + 1) (2L + 1) x/L (2L + l)jy ^)TLLM{r). (B.34) 

^B.34fl . together with the definitions of S and f, cf. Eq. (|E29}, yields 


,ijjr, 


j»(e) ■ ) ci'e« -i x/l{lTi)jl{—) fLLMir) 


X 


(/ip + (-1)^ fin) S + {fip- (-1)^ fin) t 


(B.35) 
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Including all prefactors, we get 


« OQ L 

= - V 2 >r( 2 L + {(^^+ (_i)i^„) 

^ '_LI L=1 M=—L 


k,X=±l 


X [11 ® S']lm + (/Up - (-l)^/in) [Yl ■ afc,A'^M,A - 


(B.36) 


making use of Eq. ()B.30jl once more. 

We are now ready to calculate the amplitudes and The results are given 

in Ref. [2H]; together with the amplitudes and which do not contain the 

gradient part of the photon held. Nevertheless, these contributions are strong, cf. Fig. 
due to the numerically large factor (/Xp — /x„)^ Ri 22 which appears in those amplitudes. As 
the terms with both photons coupling to the current Eq. (12. 2(1 1) . are not supported 
by this factor, these contributions are tiny and are not considered in our work. The mixed 
amplitudes also turn out negligibly small, cf. Ref. j2HI- 

So far we (explicitly) only considered one-body currents. However, there are also non- 
negligible contributions from pion-exchange currents. Fig. [Tj The corresponding expressions 
in coordinate-space representation are (with xi, X 2 the position of nucleon 1 and 2, respec¬ 
tively) 


'^s^at(C; Xi,X 2 ) = ^ (ri X ra)^ ai 5(fi - 0 ( 0^2 ■ r) + <72 5(^2 - 0 (<?i ' 


ef 


m. 


d e" 


. dr r 


for the Kroll-Ruderman (pair) current (Fig. [TKa)), and 


(B.37) 


Jl°^{^-,xuX 2 ) = (ri X 72), ■ (Vi - V2) (cfi ■ Vi) (a 2 ■ V2) - - -^^ 

\ / 777,^ — 4 | rUn \X2 — 

(B.38) 

for the so-called pion-pole current (Fig. [7Kb)), cf. e.g. The relative vector r is dehned 

as r = Ti — X 2 - 

Our numerical evaluations show that the explicit inclusion of the pole current is well 
negligible in the process and energies under consideration. Therefore we are only concerned 
with the Kroll-Ruderman current, Eq. ()B.37j) . This expression, however, is derived in the 
limit of static nucleons (denoted by the index ’stat’), i.e. the correction due to the photon 
energy is neglected. This being a rather crude approximation for uj ~ 100 MeV, which is 
close to the pion mass, we use 


fi, X 2 ) = ^ {fi X 72 ), [<71 (5(Ti - 0 (<72 ■r)+a2 <5(^2 - 0 (^i ■ r)! 


mt 


dr' 


(B.39) 


m 


^We note that our convention for the pion-nucleon coupling differs by a factor Itt from that used in 
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instead of Eq. (IB.37I1 . The function /^^(r) depends on the photon energy and is dehned as 


/f^(r) = 


2 r 


-27r 


^iq-(x2-xi) 

(27r)3 ^ {u - ^yml + ^) 


for an s-channel diagram and 


fu^ir) = 


2r 


271 


Qiq-{x2-xi) 

(27r)3 ^ (-a; - ^/ml + q'^) 


(B.40) 


(B.41) 


for the M-channel, cf. Ref. j2HI- Note that for a; = 0, Eqs. ()B.40fl and ()B.41fl reduce to the 
above expression /^^(r) = 

Now we have another current at hand, which we can use to replace J{^) in Eq. (I2.4|l . 
However, one has to be careful in order not to double-count certain contributions; e.g. 
it is not allowed to combine hf™* = —J at one vertex with = 

~ I the other, as this contribution is already included - at least partly - 

in the dominant terms, due to the use of Siegert’s theorem, cf. discussion around Eq. (ItTTI) . 
The Kroll-Ruderman current changes isospin, i.e. = — f ■ A(^) transforms 

the isospin-0 deuteron to an isospin-1 object. Therefore we need another isospin-changing 
interaction at the second vertex. Pauli’s principle guarantees that the total wave function of 
the two-nucleon system has to be antisymmetric under the exchange of the two constituents. 
Stated differently, the wave function has to fulhll (—1)'^+^+'^ = —1, i.e. in order to have T = 
1 we need S'-l-L even. The operator that turned out to be the most important one numerically 
is [Yn (5?)f]i. cf. Eq. (IB.311). The same observation has been made in uni. Nevertheless, also 
the operator Yi, which stems from cpf, gives non-negligible contributions. However, in 
the amplitudes including 0j or at the non-KR vertex, one is not allowed to use the full 
photon held in = —f ■ A(^) as explained after Eq. ()B.39|) . Therefore 

we are left with and where ’KRfull’ denotes the integral over 

the Kroll-Ruderman current, multiplied by the full photon held. There is no danger of 
double-counting as we only take into account the operator [Yl_i 0 This 

operator, however, changes the deuteron spin, whereas the matrix elements arising from 
are spin-conserving. Further contributions, like the one where J(^) = both 

vertices, turned out to be small. The evaluation of the Kroll-Ruderman terms is given in 
Ref. j2Hl5 along with further details and the analytic results for all amplitudes. 


C Total Deuteron-Photodisintegration Cross Section 

Besides complying with the low-energy theorem, cf. Section r2. 2. 41 another important check 
on our calculation is the extraction of the total deuteron-photodisintegration cross section 
from the Compton amplitude via the optical theorem. This process has been studied more 
extensively - experimentally as well as theoretically, see e.g. [33] ~ than elastic deuteron 
Compton scattering and there is plenty of data below 100 MeV to compare with. 
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The optical theorem in our normalization reads 

= - ■ s E = 0)1. (C.i) 

OJ D ^ ■' 

*=/ 

i.e. the total cross section is the sum over the imaginary part of all deuteron Compton 
amplitudes in the forward direction with identical initial and hnal photon and deuteron 
states (A/ = Aj, Mf = M*), divided by the photon energy uj. Like the elastic deuteron 
Compton cross section, this sum is divided by 6, as we have to average over the initial 
states. 

We calculate this cross section in the lab frame, in order to be able to compare to 

data. The rest of our work is performed in the yd-cm frame, which has also been chosen in 

Ref. j2Hl5 and we refer to this reference for the hnal results of our amplitudes. Fortunately, 

these amplitudes are easily transformed into the lab frame. First we note that we only 

need to sum over the s-channel diagrams, as only they become complex for photon energies 

above the deuteron binding energy B, while the w-channel amplitudes stay real for all 

nhoton energies, cf. Section 12.2.21 As the authors of Ref. uni calculate in the lab frame, 

we convert our calculation according to their work. In the s-channel the only change is 
2 2 

^ <->■ cn — because in the lab frame, the deuteron’s initial kinetic energy vanishes, 

whereas the total intermediate momentum is Pc = ki- In the cm frame we have Pi = —ki 
and Pc = 0 in the s-channel. 

Our result for the total deuteron-photodisintegration cross section is shown in Fig. 
together with data from which are described well by our calculation. In the lower 

left panel the low-energy regime is enlarged, in order to emphasize the non-vanishing value 
at threshold. The by far most important contribution at threshold stems from the singlet 
Ml-transition of which is the amplitude with the magnetic part of the photon held 

A'P, Eq. (Toil . coupling to the spin current at both vertices. It corresponds to the operator 
[Yq 0 t]i, cf. Eq. flB.31fl . which transforms the deuteron into a (singlet) Sc = 0-state. Ml 
is the shorthand notation for the magnetic coupling of a photon with L = 1. Nevertheless, 
as is well-known, already 1 MeV above threshold the cross section is completely dominated 
by the amplitude ^ (Eq. (j2.18|l i. where for L,P = 1 we have an El-interaction at 
each vertex, and this dominance holds for all higher energies. In the lower right panel of 
Fig. uni we show these two (most important) contributions to the total photodisintegration 
cross section, denoted as ’EF and ’Ml’. We observe the well-known rise of ’Ml’, as u 
approaches the breakup threshold, cf. e.g. I55j . whereas ’El’ is zero for uj = B. Note 
that El not only consists of but of all amplitudes with an El-interaction at the 

vertex of the incoming photon. 

Strictly speaking there are also contributions from the one-body current cf. 

Eq. (12.2(11). The corresnondinEr amnlitudes are ffiven in uni but are not included in this 
work, as we found that their contributions to the elastic deuteron Compton cross sections 
are tiny (of the order of 1%) and so is their effect on the total disintegration cross section. 
Nevertheless, in the high-energy regime of our calculation, say for uj ~ 100 MeV, they do 
give visible contributions to El, but these cancel nearly exactly against other terms which 
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Figure 19: Total deuteron-photodisintegration cross section derived from our deuteron 
Compton amplitudes, together with data from |1H] (open box), jlTj (diamond), jlHl (star), 
jlHl (box), jHO] (circle). The triangle corresponds to the weighted average of the data mea¬ 
sured at 2.76 MeV |51[ 1521 as determined in M- ’ El\ ’Ml’ denotes the contributions 
from the El- and the singlet Ml-transition, respectively. B is the binding energy of the 
deuteron. 


also contain Therefore, when we only look at the sum of all amplitudes contributing 
to we may well neglect the current 

We also compare our results with predictions for the strengths of electric and magnetic 
transitions close to threshold from the Effective Range Expansion given by 


_ 2 (I - 1)=^/^ 

3 7^ (|)M1 -7n)’ 


(C.2) 


mag 1 kj {l-jas + las{rs + rt)j'^-la,{r,-rt)k‘^f 

^er = 7 — - Rn) 7^-—7- , ,2 2^n -^ 

6 -\-Y [1 + k-^ a%) [1 —Tt) 

with 7 = \/rnffB. The hnal-state relative momentum is k = \pp — Pn\/2 = {uj — B), 

and for the singlet scattering length a* and the singlet (triplet) effective range {rf) we 


42 


















use tts = —23.749 fm, = 2.81 fm and rt = 1.76 fm given in j221- The explicit form of 
Eqs. (1021 ESI) is adopted from j^. 

In Fig. |2ni we compare our results with Eqs. (1021 10211 . Ending excellent agreement 
between both approaches. We also demonstrate - in the right panel - the non-negligible 
size of the KR diagrams. The only visible contributions to the total disintegration cross 
section at low energies including the Kroll-Ruderman current are the amplitudes 
which contribute to 




Figure 20: Comparison of our results (dashed) for the contributions of electric (left panel) 
and magnetic (middle and right panel) transitions to the total deuteron-photodisintegration 
cross section with predictions from the Effective Range Expansion (grey). The dotted curve 
in the right panel does not include the Kroll-Ruderman diagrams. 
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